function modelToReturn = get_model_with_added_pathway(model,pathway)

    thermodinamicsValues = evalin('base', 'thermodinamicsValues'); 
    defines = evalin('base', 'defines');
    
    [pathwayStatistics,listOfPoolMetabolites,zeroSolution, ...
    listOfPathways,listOfAllWeights, zeroSolutionBiomass, modelToReturn] = ...
        get_pathway_statistics(pathway, model, 0, true,...
        defines, thermodinamicsValues);

end